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Abstract 

The effects of the Gribov copies on the ghion and ghost propagators are investigated in 
SU (2) Euchdean Yang-Mills theory quantized in the maximal Abelian gauge. The diagonal 
component of the gluon propagator displays the characteristic Gribov type behavior. The 
off-diagonal component of the gluon propagator is found to be of the Yukawa type, with a 
dynamical mass originating from the dimension two condensate (A^A^^, which is also taken 
into account. Finally, the off-diagonal ghost propagator exhibits infrared enhancement. 
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1 Introduction 



Among the class of covariant gauges, the maximal Abelian gauge 0121 El displays several inter- 
esting features. This gauge is suitable for the study of the dual superconductivity mechanism for 
color confinement [Ij , according to which Yang- Mills theories in the low energy region should be 
described by an effective Abelian theory jSUHlEii/Sj in the presence of monopoles. A dual Meissner 
effect arising as a consequence of the condensation of these magnetic charges might give rise to 
quark confinement. Here, the Abelian configurations are identified with the diagonal components 
A^, i = 1, ...,N — 1, of the gauge field corresponding to the {N — 1) generators of the Cartan 
subgroup of SU{N). Moreover, the remaining off-diagonal components A'^, a = 1,...,N'^ — N, 
corresponding to the [N'^ — N) off-diagonal generators of SU{N), are expected to acquire a 
mass through a dynamical mechanism, thus decoupling at low energies. 

The maximal Abelian gauge can be formulated on the lattice 121131) a feature which has made 
possible to investigate the gluon propagator by numerical simulations which, in the case of 
SU{2), have reported an effective off-diagonal gluon mass of approximately 1.2GeV |niE3- 
other relevant feature of the maximal Abelian gauge is its multiplicative renormalizability to 
all orders of perturbation theory ^2 El El El- This property has allowed for a study of the 
dynamical mass generation for off-diagonal gluons, through the condensation of the operator* 
^^^^ effective potential for this operator has been evaluated in analytic form in |13j . 

providing evidence for a nonvanishing dimension two condensate 

It is worth mentioning that, although the operator has been proven to be multiplicatively 
renormalizable to all orders in the Landau, linear covariant, Curci-Ferrari and maximal Abelian 
gauges jini El El ) a satisfactory understanding of the aspects related to the gauge invariance 
of the dimension two condensate (^^) is still lacking. We refer to jl91 12U1 1211 1221 123. .24) for an 
updated analysis of this important issue. 

As other gauges, the maximal Abelian gauge is affected by the Gribov copies |2S1) whose exis- 
tence stems from a general result |27| on the lack of a globally well defined gauge fixing procedure. 
A detailed construction of an explicit example of a zero mode of the Faddeev-Popov operator 
in the maximal Abelian gauge can be found in j28j . Nevertheless, a study of the influence of 
the Gribov copies on the Green's functions of the theory in this gauge is still lacking. The 
aim of the present paper is that of providing a first analysis of the influence of the Gribov 
copies in the maximal Abelian gauge. The need for such an investigation is motivated by the 
great relevance that the Gribov copies have on the infrared behavior of Yang-Mills theories, 
as one learns from the large amount of results obtained in the Landau and Coulomb gauges 
|2nilSnil3Ilin21inSlinilESllSnil3ZlinHlinnill3l3Il. Therefore, it might be useful to improve 
as much as possible our understanding on the role of the Gribov copies in different gauges, as 

*We remind here that, due to the nonhnearity character of the maximal Abehan gauge, a slightly more 
general operator, (iy4°^Ji + ac^c"), has to be considered for renormalization purposes. The fields c", c" denote 
the off-diagonal Faddeev-Popov ghosts, while a stands for a gauge parameter. The operator (i^J^AJ^ -I- ac^c"), 
introduced in |15| . is multiplicatively renormalizable to all orders 1131 1141 125) . The maximal Abelian gauge is 
recovered in the limit a — > 0, which has to be taken after the removal of the ultraviolet divergences. Whenever 
necessary, we shall refer to 1131 for the details of the renormalization aspects of the maximal Abelian gauge as 
well as of the operator (iAJJAJ^ + ac^c"). 
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recently discussed in the case of the Hnear covariant gauges 021 • 

In the following, we shall focus on the study of the gluon and ghost propagators in the maximal 
Abelian gauge, with SU{2) as gauge group. This allows us to make a comparison with the 
results available from lattice numerical simulations. The analysis of the Gribov copies will be 
done by following Gribov's original work j26j . It turns out in fact that the construction out- 
lined by Gribov in the case of the Landau and Coulomb gauges can be essentially repeated and 
adapted to the case of the maximal Abelian gauge. We shall begin with a discussion of the gauge 
fixing condition and of the related Faddeev-Popov operator. Further, we shall generalize to the 
maximal Abelian gauge Gribov's result stating that for any field close to a horizon there is a 
gauge copy, close to the same horizon, located on the other side of the horizon^^ j26j . We shall 
proceed thus by restricting the domain of integration in the Feynman path integral to the so- 
called Gribov region, i.e. to the region in field space whose boundary is the first Gribov horizon, 
where the first vanishing eigenvalue of the Faddeev-Popov operator appears. The restriction to 
the Gribov region will be implemented by means of a no-pole condition on the ghost two-point 
function, as done in [22j. This will lead to the introduction of the Gribov parameter 7 and of 
the related gap equation, enabling us to work out the infrared behavior of the gluon and ghost 
propagators. 

A few remarks are now in order. Considering the case of the Landau gauge, it turns out 
that the restriction to the Gribov region does not eliminate all possible copies. It has been 
proven in fact that Gribov copies still exist inside the Gribov region [331 1341 I36j . To avoid 
the presence of these additional copies, a further restriction to a smaller region, known as the 
fundamental modular region, should be implemented^ . Several properties of the Gribov region 
as well as of the fundamental modular region have been established in recent years |33| I34L IHH] . 
This has been possible due to the availability of an auxiliary functional^, J-[A] = J d^x A'^A^, 
A = 1, ...,N^ — 1, whose minimization along the gauge orbit of A"^ provides a characterization 
of both Gribov and fundamental modular region. It turns out that the Gribov region can be 
defined as the set of all relative minima in field space of this auxiliary functional, while the 
fundamental modular region is identified with the set of all absolute minima of Although 
the restriction to the Gribov region does not eliminate all possible copies, its implementation in 
the Feynman path integral can be effectively worked out |3HI35). allowing one to obtain a certain 
amount of information on the infrared behavior of the gluon and ghost propagators. Such a task 
appears to be considerably difficult in the case of the modular region and, to our knowledge, it 
has not yet been accomplished. Here, a finite volume Hamiltonian approach proves to be more 
adequate [331 El HSl , see (3^1 for a review. 

Concerning now the maximal Abelian gauge, it is worth noting that a suitable auxiliary func- 
tional can be introduced also here, namely TZlA] = f d'^x A^j4^, a = 1, ...,N'^ — N, see [H 128], 
The gauge fixing condition for the off-diagonal components A^ can be obtained by requiring 
that the functional TZ[A] is stationary under gauge transformations. Moreover, a residual local 
J7(l)^~^ invariance, corresponding to the Cartan subgroup of SU{N), is still present pi I28j . 
This local invariance has to be fixed by imposing an additional condition on the diagonal com- 

^We have found useful to collect the detailed proof of this statement in Appendix A. 
"''The same conclusion holds for the Coulomb gauge. 

^The color index A runs now over all the generators of SU{N), A = 1, — 1. 
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ponents A^^ of the gauge field, which will be chosen to be of the Landau type, i.e. d^A^^ = 0. 
Analogously to the Landau and Coulomb gauges, a complete gauge fixing would require the 
implementation of the restriction of the domain of integration in the path integral to the fun- 
damental modular region for the maximal Abelian gauge, a task which is beyond our present 
capabilities. As already underlined, we shall limit ourselves to the restriction to the Gribov 
region, which turns out to correspond to field configurations which are relative minima of 7^[^]. 

The output of our results can be summarized as follows. The diagonal component of the gluon 
propagator is found to display the characteristic Gribov type behavior 

iMk)M-k)) = ^,{s,.-'-^) , (1) 

where 7 is the Gribov parameter and stands for the diagonal component of the gauge field 
in the case of SU{2), i.e. Afj, = A^. The off-diagonal propagator turns out to be of the Yukawa 
type, being given by 

(^;WAt(-.)) = (2) 

a,b = 1,2. (3) 

where m denotes the off-diagonal dynamical mass originating from the dimension two condensate 
(^A'Jj^A'Jj^ . One observes that both propagators are suppressed in the infrared. In the case of 
the ghost propagator, we find that the off-diagonal component exhibits infrared enhancement, 
namely 

Sik)\k=o ^ ^> (4) 

a 

where (c°,c°) stand for the off-diagonal Faddeev-Popov ghosts, see Appendix^ Finally, the 
diagonal component of the ghost propagator turns out to be not affected by the restriction to 
the first horizon. 



2 The gauge fixing condition for the maximal Abelian gauge 



In order to discuss the gauge fixing condition let us first remind some basic properties of the 
maximal Abelian gauge in the case of SU{2). The gauge field is decomposed into off-diagonal 
and diagonal components, according to 



An 



(5) 



where T"^, a = 1, 2, denote the off-diagonal generators of SU (2), while stands for the diagonal 
generator, 



rpa rpb 



i e'^'^T^, 
i e'^^^^ 



(6) 
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where 



Similarly, for the field strength one has 

= F;,T'' + F^,T^ , (8) 

with the off-diagonal and diagonal parts given, respectively, by 

F;, = DfAl-D-'Al, (9) 
Ffiu = d^Ai/ — duA^ + ge A^A^j , 

where the covariant derivative is defined with respect to the diagonal component A^ 

Df^d^r'-ge'^'A^. (10) 

Thus, for the Yang-Mills action in Euclidean space one obtains 

Sym = \I d^x {f;,F;, + F^^F^,) . (11) 

As it is easily checked, the classical action (|11|) is left invariant by the gauge transformations 

5AI = -D'l^'u'-ge-'Aluj, 

6A^ = -d,,u:-ge-'A1uK (12) 

The maximal Abelian gauge is obtained by demanding that the off-diagonal components of 
the gauge field obey the nonlinear condition 

<4 = ' (13) 
which follows by requiring that the auxiliary functional 

= I d^xA^Al , (14) 

is stationary with respect to the gauge transformations H12() . Moreover, as it is apparent from the 
presence of the covariant derivative D^^ , equation ()13(1 allows for a residual local U{\) invariance 
corresponding to the diagonal subgroup of SU (2) j2B]- This additional invariance has to be fixed 
by means of a suitable gauge condition on the diagonal component A^, which will be chosen to 
be of the Landau type, also adopted in lattice simulations, namely 

^^,A^ = . (15) 

Let us work out the condition for the existence of Gribov copies in the maximal Abelian gauge. In 
the case of small gauge transformations, this is easily obtained by requiring that the transformed 
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fields, eqs. (|12j) . fulfill the same gauge conditions obeyed by (^^jA^), i.e. eqs. (fTT?|) . (|T^ . Thus, 
to the first order in the gauge parameters {u!,uj^), one gets 

which, due to eas.lfT ^ .HlS p read 

M'^'^oj'' = 0, (18) 

-d^u;-ge'''d,[A';,oj') = 0, (19) 

with given by 

^ab ^ -Dl^D"^ - g'^e'^'e^'^AIA'l . (20) 

The operator Jv['^^ is recognized to be the Faddeev-Popov operator for the off-diagonal 
ghost sector, see Appendix IbI It enjoys the property of being Hermitian and, as pointed out 
in 122, is the difference of two positive semidefinite operators given, respectively, by —D'^'^D^ 
and g'^e°''^e^'^A'^A'j^. Also, one should remark that the diagonal parameter uj appears only in the 
eq. ()19|) , in a form which allows us to express it in terms of the solution of the first equation (|18() . 
More precisely, once eq. pSj) has been solved for A^, j4^, lj^, for the diagonal parameter u! one 
can write 

= -9^"'^ {K^') ■ (21) 

This feature means essentially that the diagonal parameter uj has no special role in the charac- 
terization of the Gribov copies, whose properties are encoded in ea. H18() . Also, from ea.(PT|) it 
follows that the new variable oj 

= CO + ge^'^ [a^) , (22) 

obeys 

a^cD = . (23) 

As shown in Appendix B, the change of variable ()22[) can be performed in the partition function 
expressing the Faddeev-Popov quantization of Yang- Mills theories in the maximal Abelian gauge. 
As the corresponding Jacobian turns out to be independent from the fields, transformation H22|) 
has the effect of decoupling the diagonal ghost fields from the theory. As a consequence, the 
corresponding two point function is not affected by the restriction to the Gribov region. 

3 Restriction of the domain of integration to the Gribov region 

Let us face now the implementation in the Feynman path integral of the restriction of the domain 
of integration to the Gribov region Co, defined as the set of fields fulfilling the gauge conditions 
H13|) . (|15() and for which the Faddeev-Popov operator Ai"'^ is positive definite, namely 

Co = {a„ a;, df^A, = 0, DfA'^ = 0, M^' = -D^'D'^ - g^ e'^' e''' A'^ > o} . (24) 



= 0, (16) 
= 0, (17) 
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The boundary, li, of the region Cq, where the first vanishing eigenvalue of A4°'^ appears, is called 
the first Gribov horizon. The restriction of the domain of integration to this region is supported 
by the possibility of generalizing to the maximal Abelian gauge Gribov's original result j26| 
stating that for any field located near a horizon there is a gauge copy, close to the same horizon, 
located on the other side of the horizon. We have found useful to devote the whole Appendix 1^ 
to the details of the proof of this statement. 

Thus, for the partition function of Yang-Mills theory in the maximal Abelian gauge, we write 

Z = J DA^DA^ det (m''\A)) 6 [ofAl^^ 6 (9^^^) e-^^^^V(Co) , (25) 

where the factor V{Cq) implements the restriction to the region Cq. Following j26j . the factor 
V(Co) can be accommodated for by means of a no pole condition on the off-diagonal ghost 
two-point function, given by the inverse of the Faddeev- Popov operator More precisely, 

denoting by Q{k,A) the Fourier transform of (A^"'') 



I.e. 



g(k,A) = lY.^' 



ab , ^ 



M 



ab ^ 



k) , (26) 



we shall require that Q {k, A) has no poles for a given nonvanishing value of the momentum k, 
except for a singularity at fc = 0, corresponding to the boundary of Cq, i.e. to the first Gribov 
horizon li j26(. This no pole condition can be easily understood by observing that, within the 
region Cq, the Faddeev- Popov operator A4°''^ is positive definite. This implies that its inverse, 
(TW"*) , and thus the Green function Q of eq. (j26|) . can become large only when approaching 
the horizon li, where the operator A4"^ has a zero mode. 

The Green function Q can be evaluated order by order. Repeating the same procedure of 
j26j in the case of the maximal Abelian gauge, we find that, up to the second order, 

9 V yj q q 

(27) 

where V is the Euclidean volume. We observe that the last two terms of expression ()27|) . i.e. 

Afj_(q)A^{—q) and ^Ji(<?)^^(— <?), do not depend on the external momentum k. There- 
fore, after subtraction of the corresponding ultraviolet perturbative parts^, these terms might 
yield a nonperturbative contribution to the Green function Q, corresponding to the singularity 
at A; = 0, as it is apparent from the presence of the factor l//c^ in ea. ()27() . We shall see in fact 
that these terms will give rise to a nonperturbative contribution which is proportional to the 
Gribov parameter 7. 

Thus, for Q (k, A) we shall write [201 



^At the perturbative level, these terms give rise to tadpole contributions. As such, they vanish in dimensional 
regularization, which will be implicitly employed throughout. 
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where 



^ = yT.Mq)M-9) + yJ2^:^('i)^t^^-ih (29) 

9 9 



which, in the thermodynamic hmit, V — > oo, become 
a{k,A) = g 



2k^ky f d^q A^{q)Ay{-q) 



k^ J {27rf {k-qf ' 
S = g^ j ^,A,{q)A,{-q) + ^^ j ^^Al{q)Al{-q) . (30) 

The expression for a{k,A) in ea. ()3Up can be simphfied by recalhng that, due to the Landau 
gauge condition, the Abehan component Af^{q) is transverse, namely 

q^A^iq) = . (31) 

Setting 

A^{q)A,{-q) = a;(A)(^V-^ 

uj{A) = ^Ax{q)Ax{-q) , (32) 

for a {k,A) one obtains 

„(k A) _„2feMfc.l / A Ax{q)Axi-q) f ^\ 

Note that expression (|33|) is, in practice, the same as that obtained by Gribov [2^] in the 
case of the Landau gauge. This is not surprising since a {k, A) depends only on the diagonal 
component ^^j(g), which is in fact transverse. Finally, following 26 , the no-pole condition at 
finite nonvanishing k for the Green function Q {k, A) can be stated as 

(t(0,A)<1, (34) 

with 



,N 2 1/" A Ax{q)Ax{-q) 



where use has been made of 



d\ Ax{q)Ax{-q) ( _ q^qu\ _ 3 f d\ Ax{q)Ax{-q) 
{2^f V '^^ J-^^-J (27^)' q' ' ^ ^ 

which follows from Lorentz covariance. Condition 1)34^ ensures that the Green function Q {k,A) 
in ea. (|28l) has no poles at finite nonvanishing k. The only allowed singularity is that at k = 0, 
corresponding to approaching the first Gribov horizon li. 



8 



3.1 The gluon propagator 

We are now ready to discuss the behavior of the gluon propagator when the domain of integration 
in the Feynman path integral is restricted to the region Co, eq. (|25)) . According to [22], the factor 
V(Co) implementing the restriction to Co is given by 

ViCo) = e[l-a{0,A)] , (37) 

where 6{x) stands for the step function!!. Moreover, making use of the integral representation 

/ioo+e J/- 
^gC(l-a(0,A)) ^ (38) 

for the partition function Z we get 

Z=j DA^DA,^ det (a^^^A)) exp (^C - Syu- ^ {d^'AIY - j^{d,A,f - Ca{0, A) 

(39) 

where the gauge parameters a and (3 have to be set to zero at end, i.e. a, /? — 0, to recover the 
gauge conditions ()13() . (jlSj) . In order to study the gluon propagator, it is sufficient to retain only 
the quadratic terms in expression H39|) which contribute to the two-point correlation functions 
{A''^{k)Al{-k)) and {A^{k)A,{-k)) . Thus 

^quadr = A/" ^ DA^DA^ ^^eC^-^^^ ^-^-^^-^'^(°'^)) , (40) 

where TV is a constant factor and Squadr stands for the quadratic part of the quantized Yang- Mills 
action, namely 



Therefore, recalling the expression for the factor a{0,A), eg. (|35j) . it follows 

^quadr =AA J DAf^DA^^^eKp U-l0gC-^5]^;.(g)Q^.(C,g)^.(-'Z)- ^E^;i('Z)^M-('Z)^^(-9) 



(42) 

where the quantities QfiuiC^Q) T^fiuig) are given by 



Vf^uiq) = 'Z' V - ( 1 - ^ ) 9/.^. • (43) 



9{x) = 1 for a; > 0, and e{x) = for a; < 0. 
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Note that only the factor Qp^u, corresponding to the operator appearing in the quadratic part 
for the diagonal component Ap{q) in ea. (|^^ . depends on C. Integrating over the gauge fields 
and keeping only the terms which depend on we find 

^quadr=AA J ^e'^-''^^'^ {det Q,,{Q, q))-^'\detV ^MY^ = ^' j ^/^"^^ , (44) 



where 



/(C) = C-logC-^logdet(Q^,(C,g)), 

- c-iogc-^Eiog(.^ + ^) 



(45) 



As done in expression (|44() can be now evaluated at the saddle point, namely 

^quadr ~ e^^^^) , (46) 

where Co is determined by the minimum condition 

dfiO 



dC 

which yields 



= , (47) 

C=Co 



Co 4y^g4 + M! 

Taking the thermodynamic limit, V oo, and introducing the Gribov parameter 7 [2 



2V 

we get the gap equation 



/ = Co^ ^ V^oo, (49) 



^ 2 f d Q 1 



where the term 1/Co in ea. ()48|) has been neglected in the thermodynamic limit. To obtain the 
gauge propagator, we can now go back to the expression for -Zquadr which, after substituting the 
saddle point value C = Coi becomes 

with ^ 

Evaluating the inverse of Q^uil, q) and of Vf_iu{q), and setting the gauge parameters a,P to zero, 
we get the gluon propagator for the diagonal and off diagonal components of the gauge field, 
namely 

{A,iq)AA-q)) = - ^) , (53) 
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and 

{A'^MAli-q)) = 5^^^! - Mil^ , (54) 

One sees that the diagonal component, ea. lf^ . is suppressed in the infrared, exhibiting the 
characteristic Gribov type behavior. The off-diagonal components, ea. (|54j) . remains unchanged. 
Moreover, as we shall see later, its infrared behavior turns out to be modified once the gluon 
condensate (A^^J^) is taken into account. 



3.2 The off-diagonal ghost propagator 

The off-diagonal ghost propagator can be obtained from eq.((2Hl) upon contraction of the gauge 
fields in expressions ()3U() . namely 



with 



2k^K f d^q {A^{q)A^{-q)) 



and 



B = 9^ / -0 {AM)A,{-1)} + (Aimii-q)} . (57) 

Let us consider first the factor a (k) of ea.((S^. From the expression of the diagonal propagator 
in eq. (f53|) . we obtain 

2kfj,k^ f d^q q^ ( q^q^ 



Making use of the gap equation H50() , we can write 

2k^k^ f d^q 1 / q^qy 



J (27r) q +T \ q J 

so that 

1 /,^ 2^/^^;^ / d'^q k"^ - 2kq f q^^qA 
Note that the integral in ea.()6Up is ultraviolet finite. Thus, in the infrared. A; ~ 0, one gets 



It remains now to discuss the factor B of eq. (|57)) . Making use of the dimensional regularization 
in the MS scheme, one observes that, due to the form of the off-diagonal propagator, ea. (|51)) . 
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the second term of eq.(jSZl) vanishes. Concerning now the first term, it is not difficult to see that 
it gives a contribution proportional to the Gribov parameter 7^. In fact 

{AJq)AJ-q)) = 3 / . , , , = -37 / ,4 . , = -— . (62) 



(27r)^' Mv^^ y (2^)4 ^4 + y ' y (27r)^ ^2 (^4 + y ) 32^ 

Finally, for the infrared behavior of the off-diagonal ghost propagator we have 

a(%.o-^^-^J^. (63) 
exhibiting infrared enhancement. 



4 Inclusion of the dimension two condensate (^A^A^'^ 

In this section we shall discuss the behavior of the propagators when the dimension two con- 
densate is taken into account. This condensate turns out to contribute to the gluon 
two-point function, as observed in [IHI within the operator product expansion. As such, it has 
to be taken into account when discussing the gluon propagator. 

A renormalizable effective potential for (AJ^^J^) in the maximal Abelian gauge has been con- 
structed and evaluated in analytic form in A nonvanishing condensate (j4^A^) is favoured 
since it lowers the vacuum energy. As a consequence, a dynamical tree level mass for off-diagonal 
gluons is generated. The inclusion of the condensate (^A^A^) in the present framework can be 
done along the lines outlined in j491 l5Uj . where the effects of the Gribov copies on the gluon and 
ghost propagators in the presence of the dimension two gluon condensate have been worked out 
in the Landau gauge. Let us begin by giving a brief account of the dynamical mass generation 
in the maximal Abelian gauge. Following ^j, the dynamical mass generation is accounted for 
by adding to the gauge fixed Yang-Mills action the following term 

The field a is an auxiliary field which allows one to study the condensation of the local operator 
A'^A'^. In fact, as shown in the following relation holds 

(a) = -|(yi;iA^> . (65) 

The dimensionless parameter in expression ()64() is needed to account for the ultraviolet di- 
vergences present in the vacuum correlation function {^^^(x)^^(y)). For the details of the 
renormalizability properties of the local operator ^^^^ in the maximal Abelian gauge we refer 
to [13| I14| I25| . The inclusion of the term is the starting point for evaluating the renor- 
malizable effective potential V{cr) for the auxiliary field a, obeying the renormalization group 
equations. The minimum of 1^(0") occurs for a nonvanishing vacuum expectation value of the 
auxiliary field, i.e. (cr) 7^ 0. In particular, the first order off-diagonal dynamical gluon mass 

m^ = ^, (66) 
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turns out to be |E] 

m=(^^efy Aj^^ 2.25 Ajrs- (67) 

The inclusion of the action leads to a partition function which is still plagued by the Gribov 
copies. It might be useful to note in fact that Sa is left invariant by the local gauge transforma- 
tions 

5A^ = -d^Lo-ge'^'Ay , 
5a = gA^Dfu;', 



and 



6S^ = 0. (69) 



Therefore, implementing the restriction to the region Co, for the partition function we obtain 
now 

Z = I DA^DA^ det {m''\A)) 6 (d^'^) S (d^A^) e-^'^^'+'-^Co) . (70) 

To discuss the gluon propagator we proceed as before and retain only the quadratic terms in 
expression (j7U|) which contribute to the two-point correlation functions. Expanding around the 
nonvanishing vacuum expectation value of the auxiliary field, (cr) ^ 0, one easily gets 



qua r j n ^2111 

(71) 

where the factor Q^uiC^Q) is the same as given in ea. ()43|l . while 

One sees that the inclusion of the dynamical mass m, due to the gluon condensate (^A'^j^A'^'j, 
affects only the off-diagonal sector. As a consequence, the gap equation defining the Gribov 
parameter 7, ea. (|5n|) . and the diagonal gluon propagator, ea. ((S5|l . will be not affected by the 
dynamical mass m, thus remaining the same. However, the mass m enters now the expression 
for the off-diagonal gluon propagator, which becomes of the Yukawa type, as given in expression 
p)). Note that, when the gluon condensate is taken into account, both diagonal and off-diagonal 
components of the gluon propagator are suppressed in the low momentum region. Finally, the 
infrared behavior of the ghost propagator is easily seen to display infrared enhancement 

Gik)k^o « l/k^- (73) 
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5 Comparison with lattice numerical simulations 



Having discussed the infrared behavior of the gluon and ghost propagators, as expressed by 
eqs.(^,@ and by ea. (|63|) . it is worth making a comparison with the results available from nu- 
merical lattice simulations. 



The first study of the gluon propagator on the lattice in the maximal Abelian gauge was made 
in [nj, in the case of SU{2). The gluon propagator was analysed in coordinate space and the 
Landau gauge was employed in the diagonal sector. The off-diagonal component of the gluon 
propagator was found to be short-ranged, exhibiting a Yukawa type behavior, i.e. displaying an 
exponentially suppression at large distances by an effective mass ruoff ~ 1-2 GeV. The diago- 
nal component of the gluon propagator was found to propagate over larger distances, see Fig.l 
and Fig. 2 of [Hj. These results were interpreted as evidence for the infrared Abelian dominance 
[HI lEl [Zl IS| 1 supporting the dual superconductivity picture for color confinement. 

More recently, a numerical investigation of the gluon propagator in the maximal Abelian gauge 
has been worked out in [TO]. Also here, the gauge group is SU{2) and the Landau gauge has 
been used for the diagonal sector. Moreover, the gluon propagator has been investigated now 
in momentum space, a feature which allows for a more direct comparison with our findings. 
The results obtained in [IDj show that, at low momenta, the diagonal component of the gluon 
propagator is much larger than the off-diagonal one. Several possible fits were studied for the 
components of the gluon propagator. In particular, among the two parameter fits proposed in 
[TOj, a Gribov like fit, see eq.(20) of fM, i.e. 

turns out to be suitable for the diagonal component of the gluon propagator. For off-diagonal 
gluons, a Yukawa type fit, see eq.(18) of ^Hj; ^-e. 

^o//('z) = , (75) 

seems to be well succeeded. The scalar functions, Ddiag and Doff, in eqs. lfTHl . (|7^ parametrize 
the diagonal and off-diagonal transverse components of the gluon propagator in the low momen- 
tum region 



{A^{q)A^{-q)) = Ddiagiq) 

A';,{q)At{-q)) = 5'^'Doff{q){5,.-^-i^Y (^6) 

The mass parameter rrioff appearing in the Yukawa fit ()75() is two times bigger that the corre- 
sponding mass parameter m^^ of the Gribov fit (j74|) .10) . namely 

moff ^ 2mdg , (77) 

were moff has approximately the same value as that obtained in [S], moff ~ 1.2GeV . Equation 
((77|) implies that the off-diagonal propagator is short-ranged as compared to the diagonal one. 
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Although the extrapolation of the lattice data in the region g ~ is a difficult task, which 
requires rather large lattice volumes, our results on the transverse diagonal and off-diagonal 
components of the gluon propagator can be considered in qualitative agreement with the lattice 
results, especially with the two parameter fits (|74|) and (|75j) . Concerning now the ghost propa- 
gator, to our knowledge, no lattice data are available so far. 

We remark here that the authors PI51 have also reported a nonvanishing off-diagonal longitudinal 
component of the gluon propagator which, in the low momentum region, seems to behave in a 
way similar to the off-diagonal scalar function of eq. (j75|) . Nevertheless, the analytical investiga- 
tion of this issue would require a formulation which goes beyond the original Gribov's quadratic 
approximation for the form factor a{0,A), which has been employed in the present work, see 
eas.(|35 |) . ()39() . This approximation enables us to work out a first study of the influence of the Gri- 
bov copies on the infrared behavior of the gluon and ghost propagators. Moreover, the analysis is 
by no means exhaustive and further work is certainly needed. In particular, this approximation 
does not allow to take in due account quantum corrections to the propagators in the presence 
of the Gribov horizon. One should remark in fact that the longitudinal off-diagonal propagator 
identically vanishes at the tree level, as it is easily checked from the Feynman rules stemming 
from the gauge fixing condition D'^A^^ = 0. However, due to the nonlinearity of the maximal 
Abelian gauge, one could argue that a nonvanishing off-diagonal longitudinal propagator might 
arise due to nonperturbative quantum effects. The transverse diagonal and off-diagonal propaga- 
tors, eq.(^ and eq.Q, represent a kind of first order propagators incorporating the effects of the 
Gribov horizon as well as of the dimension two condensate These propagators have to 

be used in order to investigate higher order quantum corrections as, for instance, the off-diagonal 
gluon vacuum polarization which could give rise to a longitudinal component of the off-diagonal 
propagator. Nevertheless, for a consistent evaluation of these quantum effects, we should have 
at our disposal a local and renormalizable action which takes into account the restriction to the 
Gribov region Co, ea. ()24|) . The construction of such an action has been achieved by Zwanziger 
|3H I35j in the case of the Landau gauge, where a suitable horizon function implementing the 
restriction to the Gribov horizon has been identified. Remarkably, the resulting action can be 
made local and enjoys the property of being multiplicatively renormalizable. It can be effectively 
used to evaluate quantum corrections by taking into account the restriction to the first Gribov 
horizon, see for instance the recent work [S^. Although being beyond the aim of the present 
work, we mention that the study of the horizon function for the maximal Abelian gauge is under 
investigation. Its identification would allow us to properly address the issue of the existence of 
a nonperturbative off-diagonal longitudinal gluon propagator by analytical methods. 



6 Conclusion 

In this work the effects of the Gribov copies on the gluon and ghost propagators in SU{2) Eu- 
clidean Yang-Mills theory quantized in the maximal Abelian gauge have been investigated. 

The domain of integration in the path integral has been restricted to the Gribov region Co, 
defined as the set of field configurations fulfilling the gauge conditions (|13|) . (|15)) . and for which 
the Faddeev- Popov operator M""^, ea. (|2U|) . is positive definite. Gribov's original statement 
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about closely related gauge copies located on opposite sides of a Gribov horizon has been general- 
ized to the maximal Abelian gauge, see Appendix^ providing thus a support for the restriction 
of the domain of integration to the region Cq. The dimension two gluon condensate has 
also been taken into account. 

The diagonal component of the gluon propagator displays a Gribov type behavior in the in- 
frared, eq.(^. The off-diagonal transverse component has been found to be of the Yukawa type, 
with a dynamical gluon mass originating from (^A'^A'^j^, eq.Q- Moreover, the off-diagonal ghost 
propagator exhibits infrared enhancement, eq. (j68jl . while the diagonal ghost propagator remains 
unaltered. Concerning the behavior of the transverse diagonal and off-diagonal components of 
the gluon propagator, our results can be considered in qualitative agreement with those of lattice 
numerical simulations P| llOj. 

Finally, we hope that this work will stimulate further investigation on the behavior of the 
propagators in the maximal Abelian gauge from our colleagues of the lattice community. A look 
at the off-diagonal ghost propagator would be of a certain interest for a better understanding of 
the role of the Gribov copies in this gauge. 
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A A generalization of Gribov's statement to the maximal Abelian 
gauge 

This Appendix is devoted to the generalization to the maximal Abelian gauge of Gribov's state- 
ment j20| about closely related copies located on opposite sides of a Gribov horizon. Let us 
begin by reminding that, as pointed out in |^, the Faddeev-Popov operator M"''' 

M-'{A) = -Dl'{A)D'^{A) - g^E'^''E'''AlAl , (78) 

enjoys the property of being Hermitian, being the difference of two positive semidefinite opera- 
tors given, respectively, by —D'^'^D'^ and g'^e°''^e^'^A'^^A'^. Its eigenvalues are thus real. 

Following |2n|, we can divide the space of fields fulfilling the gauge conditions (|T3|) and (|T5|) 
into regions with a definite number of bound states, i.e. negative energy solutions of the opera- 
tor M"''', see Fig.l. 

Let us look thus at the eigenvalues equation for the Faddeev-Popov operator M"'^, i.e. 

M^^i^^ = e{A)il,'' . (79) 
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Figure 1: The Gribov horizons 



For small values of the gauge fields {Af^,A'p, ea.(f7^ is solvable for positive e{A) only. More 
precisely, denoting by ei(^), e2{A), ^^{A), the eigenvalues corresponding to a given field con- 
figuration {A^,A1j^), one has that, for small all ei{A) are positive, ei{A) > 0, corre- 
sponding to field configurations for which —D^fD'j^ > e'^'^ e^'^ A'j^A'j^. However, for a sufficiently 
large value of the fields {A^, one of the eigenvalues, say ei(^), turns out to vanish, becoming 
negative as the fields increase further**. This means that the fields are large enough 
to ensure the existence of negative energy solutions, i.e. bound states. For a greater magnitude 
of (^^,^^), a second eigenvalue, say £2(^)7 will vanish, becoming negative as the fields increase 
again. Following Gribov |2n|, we may thus divide the functional space of the fields into regions 
Co, Ci, C2, Cn over which the operator M.'^^ has 0, 1, 2, n negative eigenvalues. These regions 
are separated by lines /i, ^2, In on which the operator A^"^ has zero energy solutions. The 
meaning of Fig.l is as follows. In the region Co all eigenvalues of the operator M"'^ are positive, 
i.e. M""^ > 0. At the boundary li of the region Co the first vanishing eigenvalue appears, namely 
on li the operator A4"'^ possesses a normalizable zero mode. In the region Ci the operator Ai"''' 
has one bound state, i.e. one negative energy solution. At the boundary I2, a zero eigenvalue 
reappears. In the region C2 the operator A^"'' has two bound states, i.e. two negative energy 
solutions. On ^3 a zero eigenvalue shows up again, and so on. The boundaries li, I2, I3, In, 
on which the operator Ai"'^ has zero eigenvalues are called Gribov horizons. In particular, the 
boundary li where the first vanishing eigenvalue appears is called the first horizon. See j2H] for 
an explicit example of a horizon configuration. 

**See also the argument presented in Sect. 3 of |28| . 
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It is useful to emphasize that in the region Cq, the operator Ai has only positive eigenval- 
ues. Therefore, this region can be defined as the set of all gauge fields (^^, fulfilling the 
gauge conditions eas. lfT^ . ((T3|). for which the Faddeev- Popov operator A^"* is positive definite, 
see eq.(j23]). Note also that field configurations belonging to Cq correspond to relative minima of 
the auxiliary functional 7^[^]. This follows by observing that the Faddeev- Popov operator Ai"-'' 
can be obtained by taking the second variation of TZ[A] |28j . 

Let us proceed with the generalization to the maximal Abelian gauge of Gribov's result stating 
that for any field close to a horizon there is an equivalent field, i.e. a gauge copy, located on the 
other side of the horizon, close to the same horizon, see Fig. 2 




Figure 2: The equivalent fields 



Let us start by considering a field configuration (C^, C°) located on the first Gribov horizon 
li, namely 

Df{C,)C'^ = 0, 

= 0, (80) 

where (p^ denotes a normalizable zero mode. In the following it turns out to be useful to introduce 
the diagonal component ipQ which, according to ea. ipTl) . is defined as 

^0 = -Se'^'i fc>^) . (81) 



Q2 

Let thus (^^, A^) be a field configuration located in the Gribov region Cq, close to the horizon 
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li, Fig. 2. Following j26j we write 



= C^ + a^, (82) 



where (a^, a°) have to be considered as small perturbations. The fields [A^, A'^) obey the gauge 
conditions 1)13^ and (jlSf) which, neglecting higher order terms in the small components (a^, a^), 
read 



Df{C)al-ge^'C'^a, = 0, 
5^a^ = . 

The evaluation of the energy eigenvalue e(^) of the Faddeev-Popov operator Ai°'^{A) corre- 
sponding to the field configuration (^^, A^) can be easily handled by means of perturbation 
theory, yielding 

/ d^x (2g£"-^a^Zj^^(C)(^^ + g^e'^'^e'^OjCX + ,,,, 



Proceeding as in j26j , we introduce the fields 

i/. = C/. + fl/. , (84) 



where 



a;. = - a,,<^o - se^'C^^ (85) 



have to be considered as small as compared to (C^, C^). It is not difficult to verify that, to first 
order in the small components (a^,a^), the fields (^^,A^) obey the same gauge conditions of 
namely 

D^\A)Al = 0, 

d^^A^ = 0. (86) 

The fields (A^,^°) might thus be identified with a Gribov copy of {A^,A'^), provided one is 
able to find a gauge transformation S such that 

Ay. = A^J" + A^T^ , 

Ay = A^T^^ + AyT^. (87) 
We shall look at S close to unit, in the form 

2 

S = l-a + y + 0(a3). 
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from which we obtain 



A A 1 „ \ „^o,b /ia„fe\ 9 ab„aj^bc„c , 9 /ia„a„ 

= ~ (^"^Ai" + ge A^a j - 2 " ^^^t " + ^ 



Furthermore, from eq. (|86j) . it follows 
M^\A)a^ + D';!'{A) 

-ge'^' (d,a + ge'^'A'l.a'') [o'^a' + ge'^A^a) + <7e"^4 (^|e-'^a^I)^^a^ - 



(89) 



ah Aa„b\ 9 ab„ar\bc 



9 



d„ I „_a6 /I a „ 1 ^ ^ab„a r^bc„c , ^ A<i „o,„ 
+ ge A^a ) ~ 2 "^^^j'^ 



In order to express {a, a"') in terms of Vo)' follow [inj; and set 



a" 
a 



¥'0 + 9', 

with small with respect to ((/JqiV'q). Condition H90() gives thus 

M'^\C)^' = -ge^'Dl^{C){a,^l) - ge^'^a,D^^{C)^l - g^e'^^e''\cy^ + C^)^', 



0, 

. 

(90) 

(91) 



9_^bc,^c 

2 



+ ge'^\d^^o + ge"'C';,ipt) D';{CU + ge'^C^ipo 



ge-^C 



Note that eq.(jH2) can be cast in the form 

where V"' and Q""^ are independent from ip, i.e. 

-V-{C, a, ^0) = -ge'^'D"f{C)ia,^',) - ge^^a.D'^hC)^', - g'e^'^e'\cy^ + 



(92) 
(93) 



(94) 
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and 



Equation (|98|l can be solved recursively for (f"", namely 

r = ^T''' + ^Q'''^'P' + (96) 

This allows us to obtain a recursive expression for the parameters (a, a") of eq.(jU), an thus for 
the gauge transforamtion S we are looking for, eq.(|88|). 

Moreover, recalling that 

M""^ (C) v?^ = , (97) 



one finds 



so that 



d^x ip^M^^C)^" = , (98) 



+ge^\d,ipo)D';{C)^^, - 2gh'^'e^''&^^lDf^{C)^l + g\'^' C'^C^^l^o] ■ (99) 

Let us now check on which side of the horizon li the equivalent fields (A^,^°) are located. As 
done before, we look at the energy eigenvalue e(A), given by 

^ -^^ _ J d^x {2ge--~a,D'^^{C)^l + g^£"'^e'^''(C^a^ + CX)vl) 

Finally, form eqs.ljSSI), it follows that 

e{A) = -e{A) . (101) 

Thus, if the configuration (A^,A°), close to li, is located in the region Co, e{A) > 0, there 
is an equivalent field configuration (>1^,74^), eqs. (|84|) . (|85|) . close to h, which is located in Ci, 
e{A) = — e(^) < 0. This derivation, which can be repeated to fields close to any horizon In, 
generalizes Gribov's statement to the case of the maximal Abelian gauge . 



B Faddeev-Popov quantization of Yang-Mills theory in the m2ix- 
imal Abelian gauge 

We provide here a detailed summary of the Faddeev-Popov quantization of Yang-Mills theory 
in the maximal Abelian gauge. Following let us start by giving the BRST transformations 
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of all fields, namely 



sA^ = -fd^c + ge'^'A''/ 



with 



sc" = ge'^^'c^c, 

sc = |e'*^c'^c^ 

sc" = b" , 

sc. = b , 

sb" = 0, 

sb = , (102) 
= , (103) 



where (c"^, d^) and (c, c) are the off-diagonal and diagonal Faddeev-Popov ghosts, while {b^, b) denote 
the Lagrange multipliers. For the gauge-fixed Yang-Mills theory one has 

'S'yM + 5'mAG + 5'diag , (104) 



s {Sym + 5'mAG + 'S'diag) = , (105) 

where 

Sym = d^x {F^,F^, + F^,F^,) , (106) 

and S'maGi ^diag being the gauge fixing terms corresponding to the off-diagonal and diagonal 
sectors, respectively. They are given by 

5mag = s Jd^x [ofAl + - ^ge-bri>c) 

= J d\ (b^ {DfAl + + ^DfD'^c'^ + gc'^e^' (l^j^^^^) c 

-age^^b^c^c - g^ e""^ e'"^lf A^A"^ - ^gh^'^e^^^c^c'c'^^ , (107) 

and 

5diag = s J d^xcd^A^ = j d^x (hd^A^ + cd^{d^c + ge''''Alc''y) . (108) 

Note that, for renormalizability purposes, a gauge parameter a has to be introduced in the off- 
diagonal part of the gauge fixing, ea. pU7() . The maximal Abelian gauge condition, D'j^'^A'^^ = 0, 
is recovered in the limit a — > 0, which has to be taken after the removal of the ultraviolet 
divergences |13| . In fact, some of the terms proportional to a would reappear due to radiative 
corrections, even if a = 0. See, for example, [2S1- Moreover, the action p()4|l is multiplicatively 
renormalizable to all orders of perturbation theory |121 113j . 
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Therefore, for the partition function expressing the Faddeev-Popov quantization of Yang-Mills 
theory in the maximal Abelian gauge we have 

Z = j DA^,DAlDb''DhDc''DcD-eDc e-(^YM+5MAG+Sdiag) _ (io9) 
Taking the limit a — > and integrating over the Lagrange multipliers (6", 6), one gets 

X exp I d'x ( [^DfD';y - gh-'e^WA'^A';,) + cd, [d^c + ge^'A';,c'>) ) ) . 

(110) 

To deal with the diagonal ghosts (c, c) we perform now the change of variables 

c - c = c + g^[e'^''A;c'') , 

A- ^ (111) 

all other fields remaining unchanged. It is easy to check that 

5^(^ + 5^^"'^^^')^^'^, (112) 
and that the Jacobian J corresponding to (|llip is field independent. In fact 

/ 1 ge'"'^c'' ge^^^A"^ \ 
J = det <5"^(5^i. = const. (113) 

\ 6^^^ J 

One sees thus that the transformation (|11H) allows us to decouple the diagonal ghosts from the 
theory, namely 

Z = y" DA^DAlDc^DZD-^Dc 5 (d^^A^^ 5 {d^,A^) e'^^M 

X exp [- j (fx ( (^D'^^D^^c^ - g^e^'^e^'^rc'^AlAl^ + cd"^^ ^ , (114) 

so that they can be integrated out, yielding 

Z = N j DA^DAI^Dd'Dc^ 6 {^D^A^^^ 6 (d^A^) e~^^^ 

X exp (- J d'x ( [l^DfDl^c'^ - gh-'e^'^^c'^AlAl) ) ) , (115) 

where M is an irrelevant constant factor. It is worth remarking that the change of variables 
(jlllj) seems to be a peculiar feature of the maximal Abelian gauge. One could try in fact to 
perform such a kind of transformation to decouple the Faddeev-Popov ghosts in other 
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for instance, the Landau gauge. However, it is straightforward to check now that the Jaco- 
bian of the decoupUng transformation is no more field independent. It gives back precisely the 
Faddeev-Popov determinant for the Landau gauge, so that the Faddeev-Popov ghosts show up 
again. 

Finally, integrating over the off-diagonal ghosts (c",c"), it follows 

2 = ^1 DA^DAl 6 S (d^A^) (det (^M'^') ) e"^™ , (116) 

where Ai""^ is the off-diagonal Faddeev-Popov operator, as given in eq. (j2flj) . Expression (jll6j) 
will be taken as the starting point for the implementation of the restriction to the Gribov region 
Cq for the maximal Abelian gauge. 
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